Cavity Quantum Electrodynamics of a two-level atom with modulated fields 
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We studied the interaction of a two-level atom with a frequency modulated cavity mode in an 
ideal optical cavity. The system, described by a Jay nes- Gumming Hamiltonian, gave rise to a set 
of stiff nonlinear first order equations solved numerically using implicit and semi-implicit numerical 
algorithms. We explored the evolution of the atomic system using nonlinear dynamics tools, like 
time series, phase plane, power spectral density, and Poincare sections plots, for monochromatic and 
bichromatic modulations of the cavity field. The system showed quasiperiodic and possibly chaotic 
behavior when the selected monochromatic frequencies, or ratio of the bichromatic frequencies were 
irrational (incommensurate) numbers. In addition, when the modulated frequencies were overtones 
of the Rabi frequency of the system, a single dominant frequency emerged for the system. 
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I. INTRODUCTION 

There has been a considerable amount of interest in the 
study of simple quantum mechanical systems, such as the 
simple harmonic oscillator or the two-level atom interact- 
ing with quasiperiodic electromagnetic fields. Initially, 
attention was directed to understand quantum chaos phe- 
nomenon.^^ Later on, it was observed that the atomic in- 
teraction with quasiperiodic fields could result into some 
other novel dynamical features like power broadening 
and resonance shifts,^ multiple resonance spectral struc- 
tures,^ population trapping,^ sinusoidal and square wave 
oscillation of the population^i^ and new types of dy- 
namical localization^ It was also noted that the atomic 
response to quasiperiodic fields is extremely complex, 
and the diagnostic techniques employed for characteriz- 
ing chaos could easily be utilized for studying the com- 
plexities of such interactions. 

A good number of studies are available in which inves- 
tigations of atomic systems interacting with a quasiperi- 
odic amplitude modulated field have been carried out, 
finding that the atomic density matrix elements have 
rapidly decaying autocorrelation functions. ^^ In one of 
these works, ^^ the dynamics of a two- level atom interact- 
ing with a quasiperiodic field was examined as a function 
of the field strength (which could be considered as a con- 
trolling parameter). Also, the atomic dynamics was char- 
acterized through the computation of Poincare sections 
of the Bloch vectors's motion, and the attractor's fractal 
nature. ^^ Another interesting work involved the situation 
in which a single two-level atom was interacting with a 
combination of a classical and a quantized fields in the 
cavity, so that the cavity field acted as a quantum probe 
of the dressed states defined by the classical laser field. ^^ 
Many interesting effects such as two-photon gain, cavity 
perturbed resonance fluorescence spectra, and the atomic 
squeezing in the cavity were examined. Chaotic Rabi 
oscillations under quasiperiodic perturbation were also 



studied, ^^ where a system consisting of a two- level atom 
was coupled through a time-dependent, quasiperiodic 
perturbation with two incommensurate frequencies. In 
this system, chaos was predicted by observing a rapidly 
decreasing autocorrelation and by looking at the continu- 
ous Fourier spectrum of the time-dependent physical ob- 
servable. More recently, the use of time-dependent dis- 
sipative environments to control the quantum state of 
a two- level atom was demonstrated.^^ For this case, the 
possibility to decouple dynamically the atom from its en- 
vironment leading to Markovian dynamics was studied. 

In this paper, we show an alternative approach to ex- 
plore the response of an atomic system to a quasiperi- 
odic field through its interaction with a modulated single 
mode field sustained in an ideal cavity. The frequency of 
the cavity field is modulated while its amplitude is kept 
constant, making this work different from the studies re- 
ported earlier. The interaction of a modulated field with 
a two-level atom is analyzed by studying the dynami- 
cal evolution, phase plot, Poincare section, and power 
spectrum density of the system under certain paramet- 
ric conditions. We solve the set of nonlinear differen- 
tial equations that describe the evolution of the system 
using implicit and semi-implicit algorithms specialized 
in dealing with stiff differential equations. We observe 
that for specific frequency values, the system's evolution 
is quasiperiodic, and the possibility of chaos is present. 
The cavity field is first modulated monochromatically, 
and secondly bichromatically such that the frequencies 
involved and their ratios are either irrational numbers 
or integral multiples of the Rabi frequency of the sys- 
tem. The paper is organized as follows. In section [III 
we describe the physical model of a two-level atom in- 
teracting with a modulated quantized cavity- field mode. 
Section [Till is devoted to results and discussions of the 
studied atomic system. Finally, some concluding remarks 
are presented in section [IVl 



II. MODEL AND BASIC EQUATIONS 



of three, nonlinear first order differential equations 



We consider the interaction of a two- level atom at rest, 
and atomic transition frequency cjq, with a quasiperiodic 
quantized cavity field ujc{t)- The quasiperiodic (or incom- 
mensurate) field is obtained by modulating the quantized 
cavity field frequency. The Hamiltonian of the system 
under the rotating wave approximation is 



H = ujqSz + (^c{t) ci CL -\- g{S-^a -\- a^ S- 



(1) 



where h = 1 for simplicity, a and a^ are the annihilation 
and creation operators for the cavity field, and g is the 
atom field coupling coefficient; Sz and S± are the usual 
spin-^ operators satisfying the commutation relations 



[S,,S±] = ±S±, [S+,S.] = 2S,. 



(2) 



If the modulation of the cavity field is slow, that is, the 
variation of uJc{t) with respect to ujq is small, then the 
operator a^a-\- Sz is a constant of motion, and the inter- 
action picture Hamiltonian of the system takes the form 



Hi = 5{t)Sz^g{S^a^a^S-), 



(3) 



where 6{t) = c^o — co'c(t). In this picture, the operators in- 
corporate a dependency on time, and the equation of mo- 
tion for an arbitrary time-dependent operator q is given 

by 



■iq = [q,Hi]. 



(4) 



Hence the Heisenberg's equations of motion for the 
atomic and field operators of Hamiltonian (jSj are 



ia = gS-^ 
iSz =g{S^a-a'^S-), 
iS^ = -6S^^2ga^Sz, 
iS- = SS- - 2gaSz. 



(5) 



Equation (|5]) gives rise to a nonlinear third order differ- 
ential equation in Sz 



'S z - -'S. + [5'' + Ag\N + l/2)]S, 



where N = a^a^ Sz is a constant of motion. 
If we define the following new variables 

yi = {s.), 
y2 = yi = (Sz), 
y3 = m = {Sz), 



(6) 



(7) 



then it is possible to recast (|6|) into the equivalent system 



yi 


= ^2, 




m 


= ^3, 




ys 


-^ys-is^ 


+ 4g\N+l/2)]y2 




+ V(iv 


+ i/2)^yi. 



(8) 



Solving either ([5j) or (J8]) numerically, the dynamical evo- 
lution of the system can easily be obtained for different 
parametric conditions. 

The system of nonlinear differential equations (|8j) in- 
cludes some terms {S/S and J^) that lead to rapid os- 
cillations of the solution. This makes the set of equa- 
tions stiff, so very stable numerical algorithms must be 
used to find its solutions. Because this system was small 
in size (only three differential equations) and we looked 
for solutions with moderate accuracies, we used implicit 
Runge-Kutta, and semi- implicit Rosenbrock methods to 
solve the system. ^^ 

We estimated the power spectral density (PSD) for 
each one of the studied cases to identify periodicities, 
dominant frequencies, and their correlation with the Rabi 
frequency Q of the system. To obtain these power den- 
sities, we used a simple periodogram estimator by taking 
an n-point sample of the function {Sz{t)) at equal inter- 
vals A, and using the fast Fourier transform to compute 
its forward discrete Fourier transformil 



Yk = Y.(s. 



-27rjk\/ — l/n 



A: = 0,l,2,... 



(9) 



3=0 



Then the power spectrum was defined at n/2 + 1 frequen- 
cies as 



1 



P{fk) = ^[\Yk\'^\Y^-k\'] 
fc = l,2,...,(|-l). 



(10) 



where fk took values only for zero and positive frequen- 
cies 

f' = 2^ = ^ fc = 0,l,...,-. (11) 

Finally, to facilitate the comparison and analysis of the 
dynamics of both the monochromatic and bichromatic 
modulated systems, we calculated the Poincare sections 
(PS) by taking stroboscopic views of the phase plane 
{{Sz)^ (Sz)) for just the values of 



gt 



27rn 



^/d'i+4g^N + l/2) 



n = 0,l,2,..., (12) 



corresponding to periods equal to the Rabi period of the 
systems. In this equation, Sq = S{t = 0). 



III. RESULTS AND DISCUSSION 

We solved the set of equations (|8j) numerically, using 
the implicit Runge-Kutta and semi-implicit Rosenbrock 
methods. ^^ The algorithms are adaptive (variable step- 
size /i), and the absolute and relative tolerance errors 
(atol and rtol) can be varied in accordance to the pur- 
sued accuracy, and the system of differential equations. 
Although the initial guess for the stepsize and absolute 
tolerance error we used in some of the calculations were 
relatively small {h ~ 10~^ and atol ~ 10~^), the algo- 
rithms showed stability and reliability. The results ob- 
tained for each set of parameters used were consistent 
and within the expected error. However, it is necessary 
to mention that, for monochromatic modulations, the al- 
gorithm displayed some instability and erratic behavior 
when the modulated frequencies were smaller than the 
Rabi frequency of the system. This behavior may have 
to do with the very fast oscillation of the solutions when 
the frequency differences are nearly zero, and they are 
also present as denominators in the system of differential 
equations. 

We first consider the case of no modulation of the 
quantized cavity field. In Fig. HJa), we plot {Sz{t)) as 
a function of the scaled time gt, this is the time series 
for {Sz{t)). The parameters selected were N = 3.50, 
g = 1.00, {Sz{0)) = +1/2 (the atom is initially in its ex- 
cited state), and (5o = (there is no modulation for the 
cavity field). The time series is a sinusoidal curve repre- 
senting how the population of the two- level atom evolves 
when excited with the Fock state cavity field. The math- 
ematical behavior of {Sz{t)) was quite straightforward 
in this case. Figure [TJb) shows the phase plane plot of 
{Sz{t)) against {Sz{t)). The phase trajectory is a closed 
curve with a single period of evolution. The PSD corre- 
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FIG. 1. Atom-cavity interaction for a fixed quantized cavity 
field. The parameters of the system are: N — 3.50, g — 1.00, 
^0 = 0, and (5'z(0)) = +1/2. The plots are: (a) time series, 
(b) phase plane, (c) power spectral density, and (d) Poincare 
section. 



sponding to the time series is shown in Fig. [U^c). We 
observe there is just one single frequency in the time 
series (Fig. [Il^a)) and this concur with the result ob- 
served in Fig. HJc), where we can see only one peak 
in the Fourier spectrum. This frequency matches the 
Rabi frequency of the system, which is equal to 4.0 in 
g~^ units. As expected, there was only a single point 
in the PS plot (Fig. [Hd)). Here, the PS represents 
when the phase trajectory crosses the time plane at ev- 
ery Rabi period. When we changed the parameter 5q 
to a value different from zero, all the features displayed 
in Fig. [TJa-d) remained unchanged, except that the Rabi 
frequency now turns into a more general expression given 
by [^Q +4^^(A^ + 1/2)]^/^. In addition, when we changed 
the value of N^ the generalized Rabi frequency changed. 
However, the features represented in Fig.[T]still remained 
unchanged (single frequency in the time series, closed 
curve in the phase plane, and single point in the PS). 

For the next set of parameters, we introduced mono- 
chromatic modulation of the cavity field in the form 

5 = 5q cos{ujt). 

We chose N = 2.50, g = 1.00, So = 1.00, and uj = vTf. 
The corresponding Rabi frequency was Q = 3.60 (see 
Fig. [2fa-d)). In this case, we modulated the cavity field 
with a single frequency uj = a/TT = 4.12, which is an 
irrational number. This fact gave rise to an interesting 
behavior of the dynamical evolution of the system. Fig- 
ure [21(a) shows the time series for {Sz{t)) which clearly 
exhibits a modulation pattern. The phase plot area is al- 
most filled by the phase trajectory from zero to the max- 
ima of the two variables. The PSD showed two peaks, 
one located around Q (the Rabi frequency of the system) 
and the other at uj (the modulated frequency of the cav- 
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FIG. 2. Monochromatic modulation of the cavity field. The 
time series shows a modulation pattern. Two dominant fre- 
quencies are present in the power spectrum plot. The pa- 
rameters of the system are: N = 2.50, g = 1.00, ^o = 1-00, 
(jj = vTf, and {Sz{0)) = +1/2. The plots are: (a) time series, 
(b) phase plane, (c) power spectral density, and (d) Poincare 
section. 



ity). The PS is shown in Fig. [2jd). Although there is 
a lack of synchronization between the two frequencies in 
the system, a regular star-like pattern formed, showing 
that some regions of the phase plane (states of the sys- 
tem) were visited more often than others, suggesting the 
existence of possible attractors. 

Interestingly, when uj = Q = >/T7, the situation 
changed drastically as shown in Fig. [3l We achieved this 
condition for N = 3.50. The time series now shows a sin- 
gle frequency (see Fig. [3l^a)) which is in accordance with 
the PSD of Fig. ^c). The phase plot is again a closed 
curve as depicted in Fig. [Sfb), whereas the PS shows a 
number of dots making up a wide single point in Fig.[3fd). 
These points are supposed to be at the same location; 
however, their separations are caused by the propagation 
of the computational error. In this case, when the con- 
dition cj = (^ is satisfied, the cavity modulation of the 
system synchronizes with the Rabi frequency of the sys- 
tem, and there is a single frequency or single period of 
oscillation in the dynamical evolution. We also observed 
another intriguing situation when we chose the modula- 
tion frequency to be an overtone (integral multiple) of 
the Rabi frequency, i.e., cj = ml], (m = 1, 2, 3, . . .). For 
this situation, we observed the same synchronized be- 
havior shown in Fig. [3l This phenomenon happens no 
matter the value of Q, whether it is an integer, rational, 
or irrational number. When mQ < u < {m -\- 1)Q, this 
synchronization is lost. 

In Fig. m we observed another interesting behavior in 
the dynamics of the system. This time u = {1/m) ft 
with m = 100. The ratio of cj to 1] is a fraction (with 
a positive integral denominator). Now the modulation 
frequency and the Rabi frequency are no longer synchro- 
nized with each other. Figure HJ^a), displaying the time 
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FIG. 3. The monochromatic modulation frequency of the 
cavity equals the Rabi frequency of the system. The param- 
eters of the system are: N = 3.50, g = 1.00, ^o = 1-00, 
a; = Q = Vl7, and {Sz{0)) = +1/2. The plots are: (a) 
time series, (b) phase plane, (c) power spectral density, and 
(d) Poincare section. A similar situation is presented when 
cj = ?7iO, (m — 1, 2, 3, . . .). 
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FIG. 4. Rabi frequency overtones modulating the cavity field. 
Quasiperiodicity is observed. The parameters of the system 
are: N = 2.50, g = 1.00, ^o = 1.00, u = 0.01 Q, and (S'^O)) = 
+ 1/2. The plots are: (a) time series, (b) phase plane, (c) 
power spectral density, and (d) Poincare section. 



series for {Sz{t)), apparently shows a single frequency of 
oscillation. However, the phase plot in Fig. SJ^b) reveals 
a small bandwidth around the single frequency which is 
further confirmed in the PSD in Fig.HJ^c), where there is 
no single frequency showing up in the spectrum. The PS 
for this condition is displayed in Fig. IH^d). This is typ- 
ical of quasiperiodic motion. The phase trajectory goes 
around the phase plane over and over, never exactly re- 
peating itself. As the time increases, more points will 
show up in the PS, but never coinciding with each other. 
Analyzing (J3j) for ^ = So cos(cjt), we find that the two lev- 
els in the system will cross if cos{ujt) = or t = n7r/2uj 
{n = 1,3,5,...). 

Finally, we considered bichromatic modulation of the 
cavity field by introducing the expression 

S = So[cos{ujit) + cos(cj2^)]- 

The parameters used were N = 1.50, g = 1.00, 5o = 1.00, 
uji = >/7, and UJ2 = \/T7, with a Rabi frequency Q = 3.00. 
The atom was initially in its excited state. In Fig. [5fa), 
the time series of {Sz{t)) shows what can be considered a 
sinusoidal irregular modulation. In Fig. [5fb), the phase 
plot displays quasiperiodic behavior. This is because the 
modulation frequencies selected are the square root of 
two prime numbers, 7 and 17, and the ratio of the two 
frequencies co'2/co'i, called the winding number ^^^ is irra- 
tional. Hence their combination produces incommensu- 
rate dynamical evolution of the physical variables (Sz) 
and (Sz)' The PSD corresponding to the time series is 
shown in Fig. [5fc). There are several frequencies in the 
power spectrum exhibiting the quasiperiodic behavior of 
the system. The PS plot (Fig. [5fd)) shows how the mo- 
tion of the system never exactly repeat itself, filling the 
phase plane's area as time progresses. 

Another interesting case of bichromatic modulation is 
presented in Fig. [6l For this case, the parameters used 
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FIG. 5. Bichromatic modulation of the cavity field. The time 
series plot shows an irregular pattern of modulation. Different 
frequencies determine the dynamics of the system as shown 
in the power spectrum plot. The parameters of the system 
are: N = 1.50, g = 1.00, ^o = 1-00, uJi = ^7, cj2 = ^17, 
and (*S'^(0)) = +1/2. The plots are: (a) time series, (b) phase 
plane, (c) power spectral density, and (d) Poincare section. 
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FIG. 6. Same as in Fig. [5] but now the parameters of the 
system are: N = 1.50, g = 1.00, 6o = 1.00, cji = a/TO, 
CJ2 = ^13, and {Sz{0)) = +1/2. The plots are: (a) time 
series, (b) phase plane, (c) power spectral density, and (d) 
Poincare section. 



IV. SUMMARY 



were N = 1.50, g = 1.00, 5o = 1.00, ui = ^10, and 
UJ2 = a/T3, with a Rabi frequency ft = 3.00. The time 
series is shown in Fig. [6l^a) and exhibits a similar irreg- 
ular modulation pattern to that observed in Fig. [5fa). 
Again, this is because the ratio of the modulation fre- 
quencies is an irrational number and incommensurate 
dynamical evolution emerges. The phase plot in this sit- 
uation brings out the quasiperiodic behavior of the sys- 
tem. The phase trajectory clearly outlines two distinct 
regions; one is an area densely populated by the states of 
the system, whereas the other is rarely visited by the sys- 
tem (see Fig. [6l^b)). The PSD in Fig. [6]^c) confirms this 
quasiperiodicity, where we observe a number of different 
frequencies scattered in a relatively small portion of the 
spectrum. Figure [6l^d) shows the PS of the dynamical 
system. We observe how the PS points cover the area of 
the phase plane as time moves forward. 

We also studied the behavior of the system when the 
two modulated frequencies ui and uj2 were overtones of 
the Rabi frequency Q of the system. As it happened 
before with the monochromatic case, only one dominant 
frequency emerged in the PSD plot for most of the cases. 
Two or three marginal peaks eventually appeared in the 
plots, but they were negligible compared with the single 
dominant frequency peak. We noted a slight modula- 
tion in the time series for these cases (perhaps caused by 
those marginal frequencies), but still a relative single fre- 
quency could be observed. Unfortunately, the numerical 
algorithms were not stable and reliable enough for the 
chosen parameters, so we decided not to proceed further 
in the study of this particular case. 



In summary, we presented an alternative method used 
to study the interaction of a two- level atom with a modu- 
lated quantized field of a cavity. This method comprises 
different tools used to investigate nonlinear dynamical 
systems, like the time series, phase plane, power spec- 
trum, and Poincare section plots. With these instru- 
ments, we studied the dynamical evolution of a quantum 
system. 

We described the interaction of the two- level atom with 
the cavity field using a Jaynes- Gumming Hamiltonian. 
The derived Heisenberg equations of motion of the sys- 
tem formed a set of three nonlinear first order differen- 
tial equations. Because of the stiffness and complexity of 
the equations, we used implicit Runge-Kutta, and semi- 
implicit Rosenbrock numerical methods to compute the 
solutions. A great deal of educational value is obtained 
also from this part of the work, since a simple algorithm 
for solving ODEs cannot be used here, without compro- 
mising the accuracy and reliability of the solutions. 

We explored and discussed three different scenarios: 
no modulation, monochromatic, and bichromatic modu- 
lation of the cavity field. For the different parameters 
that we used, we observed periodic and quasiperiodic be- 
havior in the population of the atomic states. We looked 
specially at those cases in which the frequencies involved 
were integer, rational, and irrational numbers, as well as 
overtones of the Rabi frequency of the system. When 
the ratio of the frequencies was an irrational number, we 
observed quasiperiodic behavior of the system. For fre- 
quencies that were overtones of the Rabi frequency of the 
system, we noted just a single dominating frequency for 
the system. 

The importance and final goal of this work are to help 
undergraduate students to understand the dynamics of a 



two-level atom interacting with a quantized mode of an 
optical cavity, when the cavity field is modulated. These 
goals are achieved by working with tools used to explore 
nonlinear dynamical systems, and solving numerically a 
set of equations of motion with specialized numerical al- 
gorithms. 
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